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The anthropic principle implies that life can emerge and be sustained only in a narrow range of 
values of fundamental constants (FCs). We extend the anthropic arguments to a regime of transient 
variations of FCs. Such regime is characteristic of clumpy dark matter models where inside the 
clumps FCs can reach values vastly different from their everyday values. We show that the passage 
of such a macroscopic clump through Earth would make Earth uninhabitable. The periodic table of 
elements is truncated and water fails to serve as a universal solvent. The anthropic principle enables 
us to improve existing astrophysical bounds on certain dark matter model couplings by five orders 


of magnitude. 


The anthropic principle implies that life as we 
know it can emerge and be sustained only in a 
certain range of fundamental constants (FCs) [1- 
3]. While the values of FCs are fixed in the Stan- 
dard Model of elementary particles, modern the- 
ories promote them to dynamic quantities. Con- 
stants are no longer constant. The variations of FCs 
have been constrained through astrophysical obser- 
vations, atomic clocks, and the Oklo natural nuclear 
reactor [4-7]. These constraints, however, presume 
slow drifts of FCs over cosmological times and are 
evaded when the variations are transient. Transient 
variations of FCs are motivated by clumpy dark mat- 
ter (DM) models such as topological defects, Q-balls 
or dark stars [8-11]. In these models, the FCs inside 
and outside the DM clumps can differ substantially. 
Encounters of Earth with DM clumps can be exceed- 
ingly rare, while the direct searches [12] so far only 
extend over a short 20-year recent history. Here, 
we show that anthropic arguments can set powerful 
constraints on FC transients over the past 4 billion 
years since the appearance of lifeforms on Earth. 

Life primarily depends on the elements in the 
first two rows of the Mendeleev periodic table. If 
a macroscopic DM clump overlaps with Earth and 
substantially affects the properties of these elements 
and the molecules containing them, it can dramati- 
cally change the conditions for life. At the quantum 
level, biology mainly depends on a set of three FCs: 
the electron mass me, the elementary charge e, and 
the Planck constant A. Variations in these FCs lead 
to a trivial isotropic scaling of the molecular geome- 
tries preserving bond angles, thus having no effects 
on biology, see Supplemental Materials (SM). How- 
ever, this changes when relativity brings in the speed 
of light c or, equivalently, the fine-structure con- 
stant a = e?/ (hc). As we show, relativistic effects 
overturn this isotropic scaling, dramatically modify- 
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ing molecular geometries. Substantial changes oc- 
cur at extreme values of a/ap = co/c ~ 10, where 
ag © 1/137 and co are the nominal values. 

Previous literature on variations of FCs predomi- 
nantly focused on small deviations from their nom- 
inal values [4, 13, 14]. Here, we explore a novel 
regime of extreme variations. We find an abundance 
of remarkable effects on the structure of atoms and 
molecules and, by extension, on the fundamental 
conditions for the emergence and sustainability of 
life. For a ~ 10ag, the periodic table shrinks to ele- 
ments from hydrogen to sulphur, the aufbau princi- 
ple of atomic structure is qualitatively changed, no- 
ble gasses are no longer inert, and water fails to serve 
as a universal solvent. If an Earth-sized DM clump 
were to have a core with œ ~ 10ag, its encounter 
with Earth would have been catastrophic. 

Qualitatively, an increase in a can be mapped into 
a reduction in the apparent value of c, amplifying 
relativistic effects. To properly describe this ultra- 
relativistic regime, we solve the four-component 
Dirac equation [15], carrying out ab initio calcula- 
tions with the DIRAC package [16] and other non- 
perturbative methods. We start with the hydrogen 
ion with a point-like nucleus. The ground state en- 
ergy is given by 


Eley, = mec? (V1— a2 -1). (1) 


As a increases, the energy is lowered towards the 
lower continuum (Dirac sea, Fig. 1) until a reaches 
1 where Elsi = —m,c?. For a > 1, the argument 
of the square root becomes negative and €15,,, aC- 
quires an imaginary part: the ground state becomes 
unstable. A more realistic analysis with a finite- 
sized nucleus shows that the 1s,/2 energy “dives” 
into the Dirac sea at —2m-c? at the critical value of 
a* = 1.04 = 143a9, see SM. For a > a*, the discrete 
1s1/2 level becomes embedded into the Dirac sea 
continuum. Then an electron-positron pair is emit- 
ted spontaneously and the vacuum becomes electri- 
cally charged, as discussed in the context of deter- 
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FIG. 1. Dependence of atomic energy levels on a. (A) 
Energy levels of atomic hydrogen as functions of a. The 
1s1/2 level dives into the Dirac sea at a* ~ 143a9. (B) 
Evolution of the aufbau principle’s sequence of shell oc- 
cupation in many-electron atoms with varying a. 


mining the critical nuclear charge with a fixed to its 
nominal value [17, 18]. For hydrogen-like ions, the 
critical value a* increases with its nuclear charge Z 
as (see SM) 


ay /a* ~ Z/168. (2) 


This remains a good approximation for multi- 
electron systems as the 1s 1/2 electrons tend to ex- 
perience the unscreened nuclear charge, with mi- 
nor correlation corrections to binding energies. In 
a molecule, a* is determined by the charge of the 
heaviest nucleus: for water and ammonia (discussed 
below) ağ ~ 18.4ao is fixed by the oxygen atom in 
H20 and aX ~% 20.9a9 by the nitrogen atom in NH3. 
Eq. (2) shows that for a given a*, only elements with 
Z S 168a0/a are stable: if a increases ten-fold, the 
entire periodic table shrinks to elements from hydro- 
gen to sulphur. 

In 1939, George Gamow published the book “Mr. 
Tompkins in Wonderland” [19] which tells a story 
about a world where FCs have values radically dif- 
ferent from those they have in the real world. In Mr. 


Tompkins’ alternative reality, where c is reduced to 
that of a speeding bicycle, c/co = ao/a ~ 4 x 1078, 
even the hydrogen atom fails to exist. 

One of the most powerful rules in atomic structure 
is the aufbau (building-up) principle which deter- 
mines the sequence in which atomic shells are filled 
with electrons. For a ~% apo, the textbook sequence 
is 181/2281/2 2p1 /22)3/2381/2 . In the ultra-relativistic 
regime (a: + a*), our calculations show that the en- 
ergy of the 3s, /2 orbital drops below that of 2p3/2, 
changing the filling order of these shells. Qualita- 
tively, the 2p3/2 orbital dives into the Dirac sea at a 
substantially larger than that for 2p; 2, see Fig. 1. 
This leads to a giant ~ Mec? fine-structure split- 
ting near the critical value a*. For the same reason, 
there is a large difference in the relativistic contrac- 
tion of the 2p3/2 and 2p;/2 shells near a*: the 2p; 2 
(and the 3s,/2) shells become submerged inside the 
2p3/2 Shell. This drives a more effective screening 
of the nuclear charge by the inner shells, causing an 
increase in the 2p3/2 orbital energies with increasing 
a in many-electron atoms, a trend opposite to that 
in the H-like ions. 

Our numerical solutions to the Dirac equation for 
second-period atoms with electrons in the 2p3/ shell 
(N, O, F, and Ne) are consistent with the described 
modified aufbau principle of Fig. 1. All the enumer- 
ated atoms demonstrate non-trivial changes in their 
ground states near their respective a*. In particular, 
the ground state of neon (Z = 10), for a S 14ap, be- 
comes 184 /2251/22P1/2381/2203/2 with the total an- 
gular momentum J = 2 instead of the nominal 
closed-shell J = 0 state. Neon is no longer inert. 
The valence-shell configuration of ultra-relativistic 
neon, 3s? /22P3 jo Closely resembles that of carbon at 


nominal a, 2s?2p?. Such neon is expected to have 
as rich a chemistry as carbon. It is intriguing to 
imagine polymers, nanostructures, and an entirely 
new biology where neon plays the traditional role of 
carbon. 

How would the extreme variation of a affect the 
structure and properties of molecules? Here we ad- 
dress this question by focusing on water and am- 
monia. All known lifeforms use water as a univer- 
sal solvent for various chemicals and as an essential 
component of many metabolic processes [20]. Am- 
monia is an important source of nitrogen required for 
the synthesis of amino acids and for building pro- 
teins in living systems [21]. At the nominal value 
of a, the water molecule is bent and the ammonia 
molecule has a tetrahedral shape. These geometries 
lead to non-zero dipole moments, which are critical 
for the chemical and biological properties of these 
molecules. At a ~ 18a 9, the water molecule be- 
comes linear and as such has a vanishing dipole mo- 
ment; water fails to serve as universal solvent. Our 


ab initio calculations show that when a increases, 
the geometry of the water molecule changes drasti- 
cally (Fig. 2 A — C). At a ~ 14a9, the bond angle 
contracts from 104.5° to 90°. At larger a but still 
below a@ ~ 18.4a9, the water molecule becomes lin- 
ear, consistent with Ref. [22]. These changes in the 
molecular geometry are induced by the relativistic 
stabilization of the 28/2 and 3s1/2 orbitals with re- 
spect to 2p; orbitals of the oxygen atom, and by the 
increased fine-structure splitting between 2p1/2 and 
2p3/2 orbitals. The formation of molecular orbitals 
(MOs) from atomic orbitals requires energy reso- 
nances and overlaps between the constituent atomic 
orbitals. Amplified relativity affects the resonances 
via the changes in energies and the overlaps via the 
varying degrees of contraction of atomic orbitals, see 
SM. 

The changes in the molecular geometry of water 
can also be understood in the context of the text- 
book valence-shell electron-pair repulsion (VSEPR) 
model: 


e At nominal a, the valence 2s and 2p atomic 
orbitals of oxygen mix to form four equiv- 
alent hybrid orbitals (sp? hybridization, or 
281 /22p1/22p3/2 in relativistic notation). Two 
of the hybrid orbitals overlap with the hydro- 
gen 1s orbitals, and the remaining two hold the 
two lone electron pairs. The repulsion between 
the four electron pairs on the hybrid orbitals 
leads to the slightly distorted tetrahedral ar- 
rangement, corresponding to the bond angle 
of 104.5°. 


e At the intermediate value a ~ 14a, the sta- 
bilization of the 2s,/2 and 2p,/2 orbitals in 
oxygen breaks down the sp? hybridization. 


The striking changes in the geometries of the wa- 
ter molecule at increased @ would lead to alterna- 
tive chemistry and biology. For example, in con- 
trast to the bent water molecules that form three- 
dimensional networks of hydrogen bonds, the ultra- 
relativistic linear water molecules could only form 
two-dimensional networks. This is anticipated to 
substantially decrease the freezing and boiling points 
of water[23, 24]. As a result, an encounter with a 
DM clump wherein a ~ a* would be catastrophic 
for life on Earth. The above observation trans- 
lates to an anthropic requirement that since life first 
emerged on its surface, Earth has not encountered a 
“large-a” DM clump. For concreteness, we assume 


This results in an energetically isolated 2p3/2 
orbital, which can accommodate up to four 
electrons forming the oxygen-hydrogen bonds. 
This stabilizes the molecular structure with 
the 90° bond angle [22]. 


e At larger a © 18a9, the lowering of the 351/2 
energy and the raising of the 2p3/2 energy 
in oxygen leads to these two orbitals becom- 
ing quasi-degenerate. This induces the ps hy- 
bridization between the half-filled 2p3;2 and 
the 3s,/2 shells [22], resulting in the linear 
molecular geometry. 


The modifications of the molecular geometry of 
ammonia at increased a resemble those of water and 
can be explained by similarity of changes in the elec- 
tronic structure of the nitrogen and oxygen atoms. 
At ao, the VSEPR model predicts the formation of 
four equivalent sp? hybrid orbitals in nitrogen, form- 
ing a tetrahedral arrangement. Three of these hy- 
brid orbitals form chemical bonds with the hydrogen 
atoms, while the fourth orbital holds the lone elec- 
tron pair, leading to the trigonal pyramidal geome- 
try, with the bond angle of 106.8°. As a increases, 
the NH3 molecule becomes first more pyramidal, but 
then completely flattens (Fig. 2 D - F). Ata = 18a9, 
the NH3 bond angle reduces to its minimum value of 
87° due to the loss of the hybridization caused by the 
relativistic stabilization of the 28,/2 atomic orbital 
of nitrogen. At a ~ 20a, the ammonia molecule be- 
comes planar with a 120° bond angle due to the in- 
duced quasi-degeneracy between the 2p3/2 and 38} /2 
orbitals, which leads to the formation of the p*s hy- 
brid orbitals. 


that the “large-a” variations are driven by a conven- 
tional quadratic coupling of photons to DM in the 
form [10, 25] 


Lom-sm > ¢°F?,/ (445) , (3) 


where ¢ is a scalar/pseudoscalar DM field, F,» is the 
Faraday tensor, and A, is an unknown characteristic 
energy scale. Inside a DM clump, the interaction (3) 
redefines a 


a= ao (1 i Prax M2)” , (4) 


where max is the field amplitude inside the clump 
(by construction, 6 > 0 outside). For a > ao, 
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FIG. 2. Molecular geometries and orbital diagrams of water (A — C) and ammonia (D — F) at different a. At 
nominal a, the orbitals are labelled and colored according to their irreducible representations in the C2, (water) and 
C3, (ammonia) point groups. At a > ao, the spin-orbital notation is used for atomic orbitals, and only the occupied 
molecular orbitals are shown. 


Eq. (4) implies A, © fmax. The time 7 between is [10] 
subsequent encounters of Earth with the DM clumps 
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FIG. 3. Comparison between anthropic and astrophysi- 
cal constraints on the energy scale Ay. 


where galactic velocity vg œ% 300 km/s, DM energy 
density ppm ~ 0.3 GeV/cm?, and d is the clump 
size. Requiring that 7 exceeds the age of life on 
Earth, we obtain 


where Dg œ~ 13000 km is Earth’s diameter. The 
anthropic constraint (6) is five orders of magnitude 
more stringent than astrophysical bounds from stel- 
lar cooling [25], see Fig. 3. 


A comparable constraint on A, comes from the re- 
quirement that nucleons remain stable. Nominally, 
proton is lighter than neutron, precluding proton de- 
cay. The mass difference results from a delicate can- 
cellation between electromagnetic (x œ) and quark 
mass splitting effects [26]. Increasing a to ~ 10ao 
spoils this balance leading both to radioactivity and 
atomic electron capture by protons. Finally, our 
predicted changes in the geometry of ammonia, one 
of the most abundant polyatomic molecules in the 
interstellar space, can be potentially used in as- 
trophysical searches for DM clumps: as a clump 
sweeps through an interstellar cloud, it would in- 
duce changes in spectral properties of the ammonia 
molecules 
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Supplemental Material 1: Invariance of 
molecular geometry under variation of 
fundamental constants in the non-relativistic 
Born-Oppenheimer approximation 


The goal of this section is to prove that the varia- 
tion of fundamental constants cause all the molecu- 
lar bonds to stretch/dilate by the very same scal- 
ing factor, leaving the angles between molecular 
bonds unaffected, see Fig. 51.4. This statement 
holds only with the assumption of (i) the non- 
relativistic approximation, (ii) infinitely-heavy nu- 
clei (Born-Oppenheimer approximation) and, (iii) 
point-like spinless nuclei. If any one of these as- 
sumptions is violated, molecular bond angles would 
vary with changing FCs. 


LS 


FIG. $1.4. Scaling all the sides of the triangle by the 
same numerical factor does not affect the value of angle 
0 (or of any other angle in the triangle). This exam- 
ple can be generalized to 3D geometry: angles and thus 
the molecular geometry are not affected by scaling of all 
the inter-nuclear distances by the same factor (isotropic 
scaling transformation). 


We begin by considering an arbitrary molecule 
containing N, point-like nuclei and Ne electrons. 
Under the enumerated assumptions, the non- 
relativistic Born-Oppenheimer (NR-BO) Hamilto- 
nian has the form 


Hon-so = -yAn +3 De 
era: z 2Me a 2 |r: = 75 
ea 


y Zae 
in [Rn E ri l 
(S1.1) 


Here, for clarity, we retained all the fundamental 
constants (FCs) and we use the Gaussian system of 
electromagnetic units. We labeled the positions of 
electrons as r; and those of nuclei as Rn, with Zn be- 
ing nuclear charges. All the terms in Hynr-pgo have 
their usual meaning: the kinetic energy of electrons, 
the electron repulsion, the nuclear repulsion, and 
the electron-nucleus attraction, respectively. The as- 
sumption of infinitely heavy nuclei allows the nuclear 
kinetic terms to be neglected. 

To determine the molecular geometry in the BO 
approximation, one first solves the time-independent 


1 Zn Zn! e? 
Ta e Ra R 


Schrödinger equation the positions of the nuclei 
fixed, 


Hyr—Bo (re|Rn) Y (re|Rn) =E (Ra) y (re[Rn) : 
(S1.2) 
Here the Hamiltonian H and thus the eigenfunc- 
tions VW and energies E depend on fundamental con- 
stants: E(R,|me_,h,e). After the potential sur- 
faces E (Ry|me, fi, e) are obtained as functions of nu- 
clear coordinates, the equilibrium nuclear positions, 
{R&4}, are determined by minimizing the energy 


min E (R,|mMe, ñ, e) 


> {R}. 
{Ra} 


(S1.3) 
We would like to explicitly factor out the de- 
pendence on FCs from Eq. (S1.2). To this end, 
we rescale all the coordinates by the same factor 
E: ri > épi, Rn > EPn. Upon substitution into 
Hynr-Bo, the kinetic energy term is transformed 
into — 3)" A2 Ap; /(2Me£?) and all the electro-static 
interaction potentials are divided by €. We may 
choose € so that the dimensionful prefactors in the 
kinetic and potential energy contributions equal one 
another, i.e., 
h2 e2 
—; =. S1.4 
meai (S1.4) 
This particular choice of € enables us to factor out 
the dependence on the FCs from the Hamiltonian 
Hynr-go. Solving Eq. (S1.4) results in 
h2 
= Å — 5 S1.5 
E= =a, (S15) 
where a is the Bohr radius. 
With this choice, the Hamiltonian admits the fac- 
torization Hnr-Bo = Enh(pe|pn), where 
En = P/M = P/E = mee* R (S1.6) 
is the atomic unit of energy (Hartree) and h(p.|py) 
is a scaled Hamiltonian which does not depend on 
FCs. The solution of the eigenvalue equation 


h(pelPn)¥? (PelPn) = E (Pn) P(PelPn) -  (S1.7) 


which is the rescaled version of Eq. (S1.2), does not 
depend on FCs either. The full energy E(R,,) is 
obtained from ¢(p,,) via 


meet 


E (Rn|me, ħ, e) = Fe E(Pn) ; 


which shows that the minimum of E(R,,|m_-,h,c) is 
attained at the point 


h2 


Mee? 


Reo = 


nr 


Pris (S1.8) 


where p*‘ is obtained by finding the equilibrium po- 
sitions 


mine(p,) +> {ef}, (819) 
{Pn} 

and is thus FC-independent. 

To reiterate, in the non-relativistic Born- 


Oppenheimer approximation, as the FCs are varied 
from their nominal values, all the equilibrium posi- 
tions are scaled by the very same factor, 

R = T Ro. (S1.10) 
Here and below all the quantities with the subscript 
0 refer to the nominal values. Equation (51.10) rep- 
resents an isotropic scaling of all the coordinates by 
the same factor, leaving all molecular bond angles 
unaffected, see Fig. S1.4. 

The fact that the isotropic scaling does not affect 
angles in a molecule of arbitrary geometry can be 
formally proven as follows. Choose {Rn}n=1...., 
be (equilibrium) position vectors of all N nuclei in 
a molecule. The angle ap between a pair of these 
vectors, Ra and Ry», is given by 


64» = cos! Fe ae 


(S1.11) 
|Ra||Ro| 


where we used the conventional definition of scalar 
products and |R,,| = y (Rn: Rn) is the length of 
the vector. Should all the position vectors be scaled 
by some factor A, R, — AR, the factors of À in 
Eq. (S1.11) cancel out. Thereby, the angles between 
molecular bonds are not affected by the isotropic 
scaling. The entire molecule undergoes isotropic 
stretching or dilation as the FCs are varied. 

In addition, as follows from our derivation, all the 
electron coordinates undergo the same isotropic scal- 
ing, 

a 


re = —Te,o- 
ao 


(S1.12) 


In particular, this means that the sizes of electronic 
clouds and atoms are scaled by the same a/apo ratio. 
Another point is that all the energies (both atomic 
and molecular) are scaled by the atomic unit of en- 
ergy 

S En 
Ero 


E (S1.13) 


These observations offer a visualization: as a 
“Mr. Tompkins” clump, introduced in the main text, 
sweeps through an atom or a molecule, all the en- 
ergy levels are gently modulated and the atom or 
molecule “breathes” in accordance with the local val- 
ues of FCs. This picture is valid in the regime of suf- 
ficiently large and slow clumps. The clumps need to 


be sufficiently large, so that there are no gradients 
of FCs across the extent of the individual atom or 
molecule. The clumps have to be sufficiently slow, so 
that the induced perturbation does not cause transi- 
tions between molecular or atomic levels. Then the 
atom or molecule follows the change in FCs adiabat- 
ically. 

It is worth emphasizing that our proof heavily re- 
lies on the possibility of factoring out all the de- 
pendence on FCs in various contributions to the 
AyrR—spo Hamiltonian. If we were to the add kinetic 
energies of the nuclei to Hyr-po, our coordinate 
scaling procedure would result in the requirement 


he e2? he he 
me E ME Mye’ 


where Mn are the nuclear masses. Generically, these 
equalities can not be satisfied simultaneously by any 
choice of the scaling parameter €. 

Our factorization procedure depends on the fact 
that the Coulomb interactions in the HNrRr-BO 
Hamiltonian exhibited power-law dependence with 
respect to distances. If the nuclei have finite size, 
the Hamiltonian no longer admits simple coordinate 
scaling. Moreover, introducing nuclear properties 
(such as finite-size charge distribution or hyperfine 
interactions with nuclear moments) into the prob- 
lem brings in another FC, m,/Agcp, where mq is 
the average mass of up and down quarks and Agcp 
is the energy scale of quantum chromo-dynamics. 

Similarly, the Dirac equation does not admit fac- 
toring out all the FCs in the Hamiltonian. In- 
deed, even in the simplest case of the hydrogen atom 
with an infinitely-heavy point-like nucleus, the Dirac 
Hamiltonian contains three terms, 


(S1.14) 


2 
e 

hp = —thea-V + bme — —. (S1.15) 
r 

Since the 4 x 4 Dirac matrices œ and £ are collec- 

tions of FC-independent c-numbers, our coordinate 

scaling procedure results in the requirement 


hc ie 
E = MeC = —. 


For arbitrary values of FCs (me, ħ, e, and c), these 
equalities cannot be simultaneously satisfied. We 
conclude that relativity must lead to the break- 
down of the isotropic scaling of atomic structure and 
molecular geometry with varying FCs. Molecular 
bond angles vary with changing FCs due to rela- 
tivistic effects. 

Since the theory of quantum electrodynamics 
(QED) is built on the quantization of relativistic 
fields, field-theoretic effects also lead to the break- 
down of the isotropic scaling with varying FCs. This 


(81.16) 


can be easily understood by examining the effects of 
vacuum polarization by the nucleus [18]. In QED, 
a nucleus is immersed into a nuclear-field-polarized 
cloud of virtual pairs of particles and anti-particles. 
Vacuum polarization leads to the replacement of the 
pure Coulomb potential —Z/r of a point-like nucleus 
by the Uehling potential. The success of our fac- 
torization procedure depends on the fact that the 
Coulomb interactions in the NR-BO Hamiltonian 
exhibits a power-law dependence with respect to dis- 
tances. The Uehling potential lacks this power-law 
dependence and, thereby, does not admit factoring 
out FCs in the resulting Hamiltonian. 

To recapitulate, the isotropic scaling of molecular 
geometry preserves angles between molecular bonds 
in the NR-BO approximation. Molecules “breathe” 
with varying FCs. 

Finally, consider a thought experiment where we 
compare lengths of two rulers of different chemi- 
cal composition. Suppose at the nominal values 
of FCs both rulers have the same lengths. As the 
FCs change, both rulers are expanding/contracting 
by the same factor in the NR-BO approximation. 
The observer would not be able to tell if FCs have 
changed. The very same argument applies to transi- 
tion frequency comparisons: in the NR-BO approxi- 
mation, all the dependence of electronic energies on 
FCs is governed by the common factor of Hartree 
energy meet /h?. Corrections to the most basic NR- 
BO approximation violate this isotropic scaling law: 
the lengths of two rulers in our though experiment 
would differ for varying FCs. Similarly, the ratios 
of transition frequencies for two different atoms or 
molecules would change with varying FCs. 


Supplemental Material 2: Critical values of a 


As discussed in the main text, the well-known an- 
alytical solution of the Dirac equation for hydrogen- 
like ions with a point-like nucleus is inadequate for 
determining the critical values a* of the electromag- 
netic fine-structure constant. In this section, we 
solve this problem using the more realistic finite- 
sized nucleus model. We also discuss the critical 
values of a for multi-electron atoms and the trunca- 
tion of periodic system at reduced speed of light. 

The critical value a*, specific to an atomic energy 
level, is determined by the requirement that the en- 
ergy £ of that level becomes equal to the Dirac sea 
threshold, € = —2mec? (see the main text). Here 
and below the rest mass energy mec? is excluded 
from £e. For some simple models of the nuclear charge 
density, this problem can be solved analytically by 
generalizing the solution for a similar problem of de- 
termining the critical nuclear charge for the fixed 


nominal value of a (see, for example, Refs. [17, 18]). 

The analytical solution can be developed for a 
spherical shell-like nuclear charge density distribu- 
tion, Psheu(r) x 6(r — R), where R is the radius of 
the nuclear charge shell. Inside the nuclear shell, 
r < R, the potential is constant V(r) = —Ze?/R 
and the solution to the Dirac equation is given by the 
free particle solution regular at the origin. Outside 
the nuclear shell, the potential is of a pure Coulomb 
character, V(r) = —Ze?/r, and the solution to the 
Dirac equation is a linear combinations of the regu- 
lar and irregular relativistic Coulomb wavefunctions. 
Setting € = —2mec? and matching the inner and 
outer solutions at r = R, we find, for the nsı/2 or- 
bitals, the following transcendental equation for a* 


ws. = 2(a*Z) cot (a*Z), (S2.1) 


where € = ,/8ZR/ao, Kiv (£) is the modified Bessel 
function of the second kind (also known as the Mac- 
donald function) with index v = 2\/(a*Z)? — 1. 
K}, (£) is the derivative of Ki, (£) with respect to £. 
Equation (S2.1) has infinitely many roots, which cor- 
respond to increasing values of the principal quan- 
tum number n. The lowest value of the root deter- 
mines a* for the 1s 1/2 orbital. 

To solve Eq. (S2.1), we need to specify the nu- 
clear charge shell radius R. We make a connection 
to the more realistic nuclear charge distributions by 
noticing that for the spherical shell distribution, the 
root-mean-square (rms) radius Ryms is identical to 
R. For proton, we take the 2018 CODATA [27] rec- 
ommended value, Rrms('H) = 0.8414(19) fm. For 
heavier elements, we use an approximation [18] R ~ 
1.6 Z'/3 fm, adequate for our semi-qualitative dis- 
cussions. From Eq. (S2.1) we find the critical value 
of a for hydrogen 151/2 to be a* ~ 1.04 or, equiva- 
lently, c* © co/143. 

One may argue that the spherical shell approx- 
imation for the nuclear charge distributions used 
in deriving Eq. (52.1) is not realistic. To address 
this question, we solved the Dirac equation for hy- 
drogen numerically using the finite-difference tech- 
niques [28]; for a uniform nuclear charge distribution 
we find a* = 1.042. This is to be compared to the 
spherical shell result of 1.040. A similar exercise 
for fermium (Z = 100, A = 257, Rims = 7.1717 fm) 
shows that the 18/2 value of a* is 1.20 x 107? for 
the spherical shell distribution and 1.18 x 107? for 
both the uniform and the Fermi nuclear charge dis- 
tributions. Thus the accuracy of the spherical shell 
approximation for the nuclear charge distribution is 
~ 0.1 — 1%, which is adequate for our goals. 

We note that the DIRAC19 package internally 
uses the Gaussian nuclear charge distributions with 


Rims given by the fitting formula [29] 
Ryms = 0.836 A13 + 0.570 (+0.05) fm. 


(S2.2) 


For a given charge Z we use the mass number A 
for the most abundant isotope. This formula re- 
sults in the proton Ryms(!H) = 1.406 fm which is al- 
most twice as large as the CODATA recommended 
value, Rrms(tH) = 0.8414(19) fm. The simple rea- 
son for this discrepancy is that Eq. (52.2) is a fit for 
atomic mass numbers A > 9, see Ref. [29]. If we use 
the value Ryms(tH) = 1.406 fm, Eq. (S2.1) results 
in a* (4H) = 1.044, slightly larger than the value of 
1.040 obtained with the CODATA Rims(tH). 

The results of our calculations for a* as a func- 
tion of Z are shown in Fig. $2.5. In this plot, the 
red curve represents the relationship between ap /a* 
and the nuclear charge Z of a finite-size nucleus with 
a Fermi charge distribution. The blue curve repre- 
sents the same dependence but for a point-like nu- 
cleus, where a* = 1/Z. This parameter space can be 
interpreted as a phase diagram: any point (ag/a, Z) 
lying above the red curve corresponds to unstable 
Dirac sea, where the 15/2 orbital is embedded into 
the Dirac sea continuum. 
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FIG. $2.5. Phase diagram of periodic system of elements 
as a function of varying fine-structure constant. Red 
curve are the results for critical values of ao/a™ as a 
function of nuclear charge Z for finite-size nuclei. Blue 
curve is the same dependence but for point-like nuclei. 
To borrow an analogy from condensed matter physics, œ 
(or c) is an order parameter that governs phase transi- 
tions. 


The finite-size nuclei critical curve exhibits a 
nearly linear dependence with a fit, 


Zmax © 168a0/Q, (S2.3) 

or, equivalently, 
ag/a* = Z/168. (S2.4) 
a < 1.23/Z (S2.5) 


The linearity of these equations can be understood 
by examining the graphical solution of the tran- 
scendental Eq. (S2.1), see Fig. 52.6 for hydrogen; 
plots for heavier elements are similar. Even with- 
out solving the equation (S2.1), it is apparent that 
the critical value of a for the 1s1/2 orbital occurs 
in the vicinity of the first zero of Macdonald func- 
tion Ki (£), where the l.h.s. approaches the ver- 
tical asymptote. The first zero of Ki (£) is given 
by Iné ~ —r/v + 1n2 — Ypuler, where YEuler = 


0.5772156649... is the Euler constant [30]. This 
leads to an analytical estimate 
1 2 1 
as Z (: p z) ; 
8 (YEuier + 4 n(V2ZR/ao)) 
(52.6) 
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FIG. $2.6. Graphical determination of critical values 
of fine-structure constant for hydrogen (Z = 1,R = 
0.8414fm). The r.h.s and the l.h.s of transcendental 
Eq. (52.1) are drawn as red and blue curves, respec- 
tively. The values of a at the intersection of two curves 
are critical values a* of a. The lowest a” is the critical 
value for 181/2, next lowest a” is for 2s1/2 and so on. 


The leading 1/Z term in Eq. (S2.6) can be recog- 
nized as the critical value for the point-like nucleus 
(see the main text). We use the smaller sign (<) be- 
cause the true value of a* lies below this asymptotic 
estimate, see Fig. 52.6. The second term is due to 
the finite size of the nucleus, with the nuclear ra- 
dius R x Z'/3. The fractional contribution of this 
second term has a weak logarithmic dependence on 
the nuclear charge, In(Z*/*), explaining the nearly 
linear dependence of the maximum allowed nuclear 
charge in Fig. $2.5. In the approximate Eq. (S2.6) 
we also restored the Bohr radius ag, showing that the 
dominant dependence is the ratio of the nuclear ra- 
dius R to the characteristic size of the atomic orbital 
ao/Z. The approximation (52.6) tends to overesti- 
mate a*. Its relative accuracy ranges from 2% for 
hydrogen to 50% for Fermium (Z = 100) as follows 


from a comparison with our numerical results. It is 
worth noting also that QED corrections to the 151/2 
energy (vacuum polarization and self-energy) tend 
to cancel [18], leaving the critical value a* largely 
unaffected by these corrections. 

For multi-electron systems, the stability of an 
atom with respect to varying FCs requires fur- 
ther qualifiers. As discussed in the main text for 
hydrogen-like ions, for a > a”, the bound 15} 2 level 
becomes embedded into the Dirac sea continuum, 
and, as such, becomes unstable, similar to Fano 
resonances in chemical physics [31]. An electron- 
positron pair is emitted spontaneously and the vac- 
uum becomes electrically charged. For low-lying en- 
ergy states of multi-electron atoms, however, the 
1s1/2 orbital is fully occupied. Then a Dirac sea elec- 
tron can not transition into the fully occupied 151/2 
orbital due to the Pauli exclusion principle [17]. Yet, 
because the rest-mass energy gap is lowered with in- 
creased a, ambient photons can promote Dirac sea 
electrons into unoccupied orbitals, i.e., Dirac sea be- 
comes unstable with respect to the interaction with 
ambient photons. 

As for the critical values of multi-electron sys- 
tems, we computed the Dirac-Hartree-Fock (DHF) 
energies of the 1s 1/2 orbitals in several atoms as a 
function of a. We find that the hydrogen-like ion 
result (S2.3) for a* remains an excellent approxi- 
mation. Indeed, the energies of the deeply-bound 
1s,/2 orbitals in atoms and molecules are strongly 
dominated by the interaction with the nucleus with 
a small corrections from the interaction with other 
electrons. 


Supplemental Material 3: Results for 
many-electron atoms 


In this section, we discuss the effects of increasing 
a on the spectrum of a many-electron atom. We 
will focus on neon. Results for other second-row 
elements may be found in Ref. [32]. 

The atomic spectrum of Ne was calculated at dif- 
ferent values of a using the Kramer Restricted Con- 
figuration Interaction (KR-CI) method, as imple- 
mented in the DIRAC19 package [33]. The Dirac- 
Hartree-Fock (DHF) calculation was first carried 
out to obtain the reference wavefunction for the 
CI step. In the DHF calculation, the 3s1/2 orbital 
was included in the active space via the average-of- 
configuration open-shell framework [34]. This en- 
sures the balanced description of atomic states at 
the CI step and allows us to assess the effects of the 
relativistic stabilization of the 3s1/2 orbital. 

At the nominal ag, the atomic states are labeled 
in the conventional L-S (Russell-Saunders) coupling 


scheme: ?S+tIL}. At increased a, however, ampli- 
fied relativity leads to the breakdown of the LS cou- 
pling scheme, as only the total angular momentum 
J = L + S is conserved [28]. Thereby, we label the 
states as J”, where J is the value of the total angu- 
lar momentum and ~ is the parity of the state. If 
there are multiple states of the same J” symmetry, 
we distinguish them by appending their sequential 
number n: J7(n), where the states are enumerated 
in the order of increasing energy. In our notation 
for electronic configurations, for brevity, we suppress 
the 1s,/2 and 2s1/2 shells, as these remain always 
doubly occupied for our considered low-lying energy 
states. 

Since the typical distance of an electron from the 
nucleus decreases with increasing a due to relativis- 
tic contraction, the basis set used in our calculation 
needed to be calibrated to accurately describe the 
electronic density near the nucleus at increased a. 
The calibrating procedure was carried out by con- 
sidering the hydrogen-like ions of Ne as follows. We 
gradually increased a until the 1s, 2 ground state of 
Ne? dived into the Dirac sea. The size of the basis 
set and the largest exponents were chosen to match 
the critical values of a* obtained using such a basis 
set with that predicted by solving the transcendental 
Eq. (S2.1). Additionally, the validity of the basis set 
was verified by comparing the energy level orderings 
they generated with those predicted by the finite- 
difference solution of the Dirac equation. To obtain 
the correct energy ordering, it required to augment 
standard basis sets with additional p basis functions. 

For Ne, the eleven p basis functions in the origi- 
nal unc-aug-cc-pV6Z basis set were augmented to a 
total of nineteen. The exponents of the new func- 
tions were obtained by subsequently multiplying the 
largest p exponent by 3. For simplicity, the same 
exponents were used for the s basis functions. This 
procedure was carried out until a match with the so- 
lution to the transcendental Eq. (S2.1) was obtained, 
while maintaining the correct ordering of the energy 
levels. The resulting modified basis thus included 
nineteen s and nineteen p basis functions with the 
largest exponent being 1.4 x 10°. Basis functions of 
higher angular momenta were left unchanged. With 
such a modified basis set, the critical a* value ob- 
tained using the DIRAC19 program matched that 
predicted by Eq. (52.1), namely, aX, ~ 14.8a9. 

Note that to avoid the collapse of many-electron 
wavefunctions into the Dirac sea [35], the so-called 
“no-pair” Hamiltonian [36, 37] was used in fully rel- 
ativistic electronic structure calculations 


a ok 1 
HT no —pait = 5 hp(i) + z D ` (53.1) 
i iAj 33 


Here the first term is a sum of the Dirac Hamilto- 
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nian hp(i) describing an electron i moving in the po- 
tential of a finite-size nucleus and the second term 
describes the Coulomb repulsion between the elec- 
trons. Notice that the e — e interaction is sand- 
wiched between projection operators A,+ which ex- 
clude states from the Dirac sea continuum of hp. 


We have described the methods with which we 
computed the low-lying energy states of Ne at dif- 
ferent values of a. Below, we present the results of 
our calculation. We find that as a is increased, the 
energies of the excited states of Ne exhibit various 
interesting features. Relative to the energy of the 
“nominal” ground state (the ground state at nom- 
inal a), an excited state energy may rise or fall in 
the regime a ~ 10ag, leading to several crossings of 
levels. However, as a is increased further, all excited 
states eventually stabilize with respect to the nomi- 
nal ground state. Even more remarkably, as œ nears 
a*, the energies of some excited states decrease so 
much that these states become the “new” ground 
states themselves. 


We show below that the electron-configuration 
picture is sufficient for the qualitatively explanation 
of these phenomena. The lynch-pins of this exposi- 
tion are the facts that as a increases, (i) the 2p1/2 
energy decreases, (ii) the 2p3/2 energy increases, (iii) 
the 3p 1/2 energy falls, but with a slower rate than 
that of 2p1/2, (iv) the 3p3/2 energy rises but with a 
slower rate than that of 2p3/2, and (v) all s1/2 en- 
ergies fall. It is worth stressing that the rise in the 
energy of 2p3/2 (and similarly of all other p3,2) or- 
bitals is only present in multi-electron atoms. In a 
hydrogen-like atom, although the fine-structure con- 
tribution to a 2p3/2 energy increases with increasing 
a, the gross-structure contribution decreases (it is 
a negative quantity whose magnitude gets larger), 
leading to an overall decline of the 2p3/2 energy. 
In a multi-electron atom, however, the contractions 
of the inner 1s1/2, 281/2, and 2pı/2 orbitals with 
increasing œ leads to more effective screening of 
the nuclear charge, thus reducing the magnitude of 
the gross-structure contribution to the 2p3/2 energy. 
This in turn means that the 2p3/2 energy increases 
with increasing a. 


With this general picture in mind, let us discuss in 
details the changes in Ne atomic spectrum as a in- 
creases. At the nominal a = ag, the closed-shell 2p°® 
ground electron configuration 2p of neon results in a 
single atomic state with even parity, ‘Sp. The 2p°3s 
first excited configuration gives rise to four excited 
states with odd parity: 7[3/2]$, ?[3/2]9, 7[1/2]8, and 
?(1/2]9. The second excited configuration, 2p°3p, 
generates a manifold of atomic states with even par- 
ity, among which we consider the three lowest states, 
?11/2]1, 7[5/2]3, and ?[5/2]2. We therefore include a 
total of eight atomic states of neon in our discus- 
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FIG. $3.7. Energy spectrum of neon atom as a function 
of a. The nominal ground state 0°('So) is used as a 
reference (blue horizontal line). At the nominal a, the 
levels are labeled using the conventional L-S coupling 
scheme [38]. At larger values of a, the levels are labeled 
as J", where J is the total angular momentum and 7 is 
the parity of the state. There is a substantial reshuffling 
of sequence of energies near the critical value of a ~ 
a* = 14.8a9. In this regime, the ground state of neon 
atom becomes the open-shell 2° state. 


sion, see Fig. $3.7 in the main text. At larger a, 
these states are labeled, in the order they are intro- 
duced above, as 0°, 2°, 1°(1), 0°, 1°(2), 1°, 3°, and 
2°, respectively. 

Throughout this discussion, the nominal ground 
state 0°, which corresponds to tSo at a = ag, is used 
as a reference; and all state energies are reported as 
the ratio (E—E¢)/|Eoe|. This choice does not imply, 
however, that the state 0° remains unaffected as œ 
increases. In the range ag < a Ê 13a9, the 0° state 
retains 0° its closed-shell configuration of 2p? /22P3 jx 
However, at a ~ 13ao, the 3s1/2 and 2p3/2 ener- 
gies become close enough so that 2p? 12203 J2351 J2 
emerges as an appreciable contribution in the CI ex- 
pansion of 0°. Note that two electrons are trans- 
ferred from 2p3/2 to 3s1/2 to preserve total angular 
momentum and parity. The appearance of the open- 
shell configuration 24 192P3 19384 9 in the CI expan- 
sion of the nominal ground state indicates that at 
a = 14a9, neon is no longer chemically inert. 

The fact that increasing a has the effect of “acti- 
vating” the naturally inert neon may also be under- 
stood by considering the state 2°, which corresponds 
to ?[3/2]§ at a = ao. The main configuration for 2° 
in the regime ao < a S 18a is 274 192P3 9381/2 
containing the orbital 3s1/2 which has the effect of 
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destabilizing 2°. As discussed in the previous para- 
graph, in the vicinity of a ~ 13a9, the ground state 
0° acquires the component 2p? J22P3 7238} jz Since 
the 3512 orbital still lies above 2p3/ in this regime, 
the 2° energy falls below that of 0° and 2° briefly 
becomes the “new” ground state of neon. However, 
as a nears 14a9, the 2p3/2 and 3s1/2 orbitals cross, 
raising 2° back above 0°. 

The state 1°(1) (nominally ?[3/2]?) displays a de- 
pendence on varying a similar to that of 2°. In 
the regime ag < œa S 13a, its CI expansion is 
dominated by 2p /22P3 9381/2 and 211/223 9381/2 
which stabilize its energy relative to the ground 
state. However, due to the second configuration 
where 2p3/2 is quadruply occupied, the energy of 
1°(1) does not lower as dramatically as that of 2°. 
In particular, the 1°(1) energy never falls below the 
ground state energy. As aœ approaches then passes 
14a9, 381/2 crosses below 2p3/2 and 1°(1) destabi- 
lizes relative to 0°. The rate of destabilization is 
reduced as œ approaches a* due to the appearance 
of the configuration 24 /92p3/2384 9481/2 in the ex- 
pansion of 1°(1). 

The last odd-parity state included in our discus- 
sion is 0° (nominally ?[1/2]$) which, for ao < a < 
10a9, comprises mainly of 291 /22P3 9381/2 which, 
similarly to the case of 1°(2), causes 0° to destabilize 
relative to the nominal ground state. However, un- 
like 1°(2), the absence of the configuration 1°(2) in 
the expansion of 0° means that its energy continues 
to rise until a ~ 10a09, where an avoided crossing 
with higher (nominally) ?[1/2]§ state (not shown) 
changes the configuration of 0° to 2% (923 /23d5/2- 
Since the 3ds5/2 orbital destabilizes with increasing 
a at a much slower rate than 2p3/2, the 0° energy 
experiences a steep downturn relative to the ground 
state energy for a = 10ap. 

Next, we consider the even-parity excited states 
1° and 2°, which correspond to 7[1/2], and ?[5/2]2 
at a = ag. These states both start out at as combi- 
nations of 274 192P3 /23P1/2 and 24 j22P3/23P3/2 As 
a increases, the 3p, /2 energy decreases and the 3p3/2 
increases at a much slower rate than that of 2p1/2. 
As a result, the states 1° and 2° generally desta- 
bilize with respect to the nominal ground state 0°. 
However, the stabilization pattern of 2° display a 
peculiar feature. At a ~ 13a0, the configuration 
2p /22P3 /23P3/2 is replaced by 24 192P3 19384 9 from 
the CI expansion of 1° and 2°. As a result, the de- 
crease of the 2° energy becomes even more precipi- 
tous and at œ & 14a, it becomes less than the 0° 
energy, thus making 2° the new ground state of neon. 

Finally, we consider the state 3°, which is labeled 
?15/2|3 at ao. The dominant configuration in this 
state remains 2% 19223 /93P3/2 for all ag < a < a*. 


As a result, it destabilizes continuously relative to 
the ground state, albeit not rapidly enough for it 
to cross the nominal ground state anywhere in the 
interval ag <a Sa’. 


Supplemental Material 4: Results for molecules 


In the relativistic picture, the states of atoms and 
molecules are described by four-component Dirac 
spinors Y = (We, Y, yg, We), where L and S indi- 
cate the large and small components, respectively, 
while œ and 8 describe the spin degrees of free- 
dom. The spinor components are, in general, com- 
plex numbers so a general collection of four such 
components has eight degrees of freedom. However, 
since the spatial and spin degrees of freedom are 
coupled, the symmetry of the Dirac spinors is de- 
scribed by the double groups, where the total spinor 
transforms under fermion irreducible representations 
spanned by half-integer spin functions [34]. Further- 
more, the real and imaginary parts of each spinor 
component are spanned by boson irreducible repre- 
sentations, which are irreducible representations of 
conventional single point groups. Therefore, each 
spinor component can be described by scalar func- 
tions, or orbitals. 

The symmetries of the water and ammonia 
molecules are described by the Co, and C3, dou- 
ble groups. For example, exploiting the symmetry 
of the Dirac Hamiltonian, it can be shown that in 
the C2, double group, the real and imaginary parts 
of the large component transform under the (a1,a2) 
and (b1,b2) boson irreducible representations for Yg 
and we, correspondingly [16]. At nominal a, molec- 
ular orbitals of ammonia and water are spanned 
only by a real or imaginary part of a single compo- 
nent, neglecting vanishing contribution from other 
components. Therefore, these orbitals are described 
by a single irreducible representation, in compli- 
ance with results from non-relativistic calculations. 
At increased a, however, the molecular orbitals are 
spanned by multiple real and imaginary parts of the 
spinor components and one cane no longer assign 
a single irreducible representation to molecular or- 
bitals. For this reason, the symmetry labels in the 
MO diagrams of ammonia and water are presented 
only at the nominal speed of light, but not for larger 
q. 
In the MO diagram of water (Fig. 2 A-C in the 
main text), the o and o* linear combinations of 1s 
orbitals of two hydrogen atoms have the a; and b2 
symmetries. For oxygen at nominal a, the atomic 
orbitals 2s;/2 and 31/2 have the symmetry of a; 
whereas the 2p orbitals have the symmetries of a1, 
bg, and bı. In the ammonia MO diagram (Fig. 2 D-F 
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in the main text), the three linear combinations p1, 
we, and w3 of 1s orbitals belonging to the three hy- 
drogen atoms have a; and doubly degenerate e sym- 
metries. At nominal a, the 2s;/2 and 3581/2 atomic 
orbitals of nitrogen have symmetry a; whereas the 
2p orbitals have symmetries a, and e. In atomic cal- 
culations of oxygen and nitrogen, the 3s1/2 spinor 
was included in the average-of-configuration Dirac- 
Hartree-Fock method to assess the effect of stabi- 
lizing higher lying spinors on molecular bonding. 
For clearer comparison of diagrams, we keep the en- 
ergy unit constant and equal to that at nominal a. 
To better demonstrate the changes in the electronic 
structure of molecules at increased a, we calculated 
the radial density distribution in DIRAC19 [33, 34] 
for each molecular orbital of water and ammonia as 


p(r) =| T p(r)r? sin 6 dô do, (S4.1) 


where p(r) = |W(r)|? is the electron density and V(r) 
is the MO wave function. 
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FIG. $4.8. Fractional changes in the bond distance r and 
the bond angle @ in the water molecule at different values 
of a. The critical value of a for oxygen is ad = 18.4a0. 


The Walsh correlation diagrams of the MO the- 
ory provide more insights into the relation between 
the electronic structure and geometry of the water 
molecule, see Fig. $4.10. These diagrams show the 
energies of valence MOs as functions of the bond 
angle. Because the total energy of a molecule can 
be approximated as the sum of MO energies, the 
Walsh diagrams can be used to predict the values of 
the bond angle that minimize the total energy. At 
nominal a, the interplay between the HOMO-1 and 
HOMO-2 energies minimizes the total energy at the 
bond angle of 104.5°. As a increases, the 21/2 or- 
bital stabilizes and its contribution to HOMO is di- 
minished, leading to the minima of the HOMO and 
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FIG. $4.9. Radial density distributions for four occu- 
pied molecular orbitals of water at different values of a. 
The oxygen atom is placed at the coordinate origin. The 
molecular orbitals are labeled in the order of increasing 
energy. The dark blue solid curve shows the distribu- 
tion for the lowest occupied molecular orbital, HOMO-3 
(HOMO stands for the highest occupied molecular or- 
bital); the dark red dashed curve - HOMO-2; the dark 
green dotted curve - HOMO-1; and the magenta dot- 
dashed curve - HOMO. For better resolution, the dis- 
tance from the oxygen atom is given in the units of the 
unscaled nominal Bohr radius. 


total energies at the 90° bond angle. At larger a, 
the relativistically stabilized 3s,/2 orbital starts to 
contribute to HOMO, leading to the linear geometry. 

The Walsh correlation diagrams for ammonia are 
shown in Fig. 54.13. At the nominal and intermedi- 
ate a, the minima of the sum of orbital energies cor- 
respond to the trigonal pyramidal molecular geome- 
tries. For even larger a, the HOMO stabilization at 
large angles dominates the energy changes of other 
molecular orbitals, leading to the planar molecular 
geometry. 


Supplemental Material 5: Dark matter clumps 
with large variations of fundamental constants 


In this section, we discuss in more details the prop- 
erties of a dark matter (DM) clump whose interior 
facilitates a large variation of the fundamental con- 
stants, in particular the fine-structure constant a. 

Depending on its mass and whether it is self- 
interacting, DM can exist in a variety of configu- 
rations. For example, DM in the ~ GeV-TeV mass 
range is often searched for in the form of weakly- 
interacting massive particles (WIMPs) in particle 
detectors. In contrast, ultralight DM behaves more 
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FIG. $4.10. Walsh correlation diagrams of water show- 
ing the energies of four valence molecular orbitals as 
functions of the bond angle at different a. In panels B 
and C, the lowest occupied molecular orbital is omitted 
as it reduces to the 2s,/2 atomic orbital of oxygen and 
does not participate in the formation of oxygen-hydrogen 
bonds. The arrows indicate the bond angle at which the 
molecular geometry stabilizes. The energies of molecu- 
lar orbitals at the equilibrium bond angles correspond to 
those in Fig. 2 A-C in the main text. 
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FIG. $4.11. Fractional changes in the bond distance r 
and the bond angle 0 in the ammonia molecule at differ- 
ent a. The critical value of a for nitrogen is ax = 20.9a0. 


like classical fields due to macroscopic occupation 
numbers. Further, if the ultralight DM field is 
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FIG. $4.12. Radial density distributions for four occu- 
pied molecular orbitals of ammonia at different a. The 
nitrogen atom is placed at the coordinate origin. The 
molecular orbitals are labeled in the order of increasing 
energy. The dark blue solid curve shows the distribu- 
tion for the lowest occupied molecular orbital, HOMO-3 
(HOMO stands for the highest occupied molecular or- 
bital); the dark red dashed curve - HOMO-2; the dark 
green dotted curve - HOMO-1; and the magenta dot- 
dashed curve - HOMO. For better resolution, the dis- 
tance from the nitrogen atom is given in the units of the 
unscaled nominal Bohr radius. 


non-self-interacting, as in the cases of the pseudo- 
scalar axions or scalar dilatons/moduli, it forms a 
uniform background DM field oscillating coherently 
at the Compton frequency (“wavy” DM). In this 
paper, we are concerned with ultralight DM fields 
which possess self-interactions and, as such, can 
form localized structures (“clumpy” DM). Examples 
of “clumpy” DM models include Q-balls [9, 39-41], 
Bose stars [42—44], topological defects [8, 45, 46], ax- 
ion quark nuggets [47—49] and “dark blobs” [50]. 

For our purpose, the specific details of the DM 
model are unimportant and we consider a generic 
ultralight DM field ¢ that forms a “gas” of clumps. 
For topological defects, the clump’s transverse size 
is of the order of the field Compton wavelength, 
d ~ h/(mgc) where mg is the mass of DM candidate 
field. An Earth-sized clump would therefore corre- 
sponds to mg ~ 10714 eV. In other models, e.g., Q- 
balls, the clump size is determined by the specifics of 
the self-interaction potential [39]. The exact relation 
between d and mg is irrelevant for our estimates, but 
the underlying candidate DM fields are necessarily 
ultralight. 

Assuming that the DM clumps saturates the lo- 
cal DM density ppm ~ 0.3 GeV/cm?, and that their 
velocities in the Solar system are vg ~ 10~%c ~ 300 
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FIG. $4.13. Walsh correlation diagrams of ammonia 
showing the energies of four valence molecular orbitals 
as functions of bond angle at different a. In panels B 
and C, the lowest occupied molecular orbital is omit- 
ted as it reduces to 2s1/2 core spin-orbital of nitrogen 
and does not participate in the formation of nitrogen- 
hydrogen bonds. The energies of molecular orbitals at 
the equilibrium bond angles correspond to those in Fig. 2 
D-F in the main text. 


km/s (~ Solar system velocity relative to the DM 
halo), the average time between successive “close en- 
counters” with the clumps is given by [10] 


max 


So ee 


55.1 
Ug ppmd he , ( ) 


where @max is the DM field amplitude inside the 
clump (¢ > 0 outside). 

Due to the coupling of DM to baryonic matter, an 
encounter with such a DM clump may cause tran- 
sient variations of fundamental constants [10]. For 
Earth-sized DM clumps with 7 of the order of a 
few years, these transient signals may be detected 
by a network of sensors searching for perturbation 
patterns as the clumps sweep through the network. 
Such direct searches are being carried out with net- 
works of atomic clocks, such as atomic clocks aboard 
satellites of the Global Positioning System [12] or the 
trans-European network of laboratory clocks [51]. 
However, if 7 reaches the cosmological scale of a 
few billion years, constraints on the DM parameters 
can be derived via anthropic arguments presented in 
the main text. In some sense, the lifeforms serve as 
DM detectors. 


In particular, invoking the anthropic argument 
that a DM clump with a fine-structure constant @ 
much larger than the nominal value ag ~ 1/137 
has not encountered Earth since the beginning of 
life on Earth is equivalent to requiring T = 4 Gy. 
Eq. (S5.1) then translates into 

Pmax Z 5 x 10° TeV. (S5.2) 

The assumed coupling of photons to DM has the 
form (see main text) 

Lpm-sm > Fi /(44°7), (55.3) 
where F,» is the Faraday tensor and A, is the char- 


acteristic energy scale. This coupling re-scales a in 
the core of the clump, 


Omax = 09 (1 — Pig / A2) . (S5.4) 
For Qmax >> Qo, the regime of our interest, Eq. (S5.4) 
implies that A, © max which, when combined with 
Eq. (S5.2), gives the limit 
A, Z 5 x 10° TeV. (S5.5) 

We note that for A, satisfying Eq. (55.5), the DM 
clumps are not subject to the screening by Earth’s 
atmosphere [52] implied by the interaction (55.3). 
The limit in Eq. (S5.5) is largely insensitive to the 
specific value of Qmax, resulting to the insensitivity 
of Eq. (S5.5) to screening effects, e.g., by the ocean 
water. 

Let us now consider the implications of Eq. (55.2). 
Assuming an Earth-sized clump, d = 13,000km, the 
energy density inside the clump is found to be 


Pinside ~ Prrax/ d © 10° kg/cm? , (S5.6) 
and its total mass-energy is 
E~ Pinside X d? 
x 3 x 10!" kg ~ 107° mass of Earth. (95.7) 
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We observe that an Earth-sized DM clump with 
max ~ 5X 10° TeV only weights as much as a small 
asteroid. In particular, such a clump does not col- 
lapse into a black hole. Further, if such a clump en- 
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that the DM clumps saturate the local DM density, 
we may estimate their number density n as 
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which is the same order as the size of the Solar sys- 
tem. 
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of an electron-positron pair when the 1s,/2 atomic 
level dives into Dirac sea. This characteristic 
sweep time is also much longer than the charac- 
teristic picosecond timescale for nuclear motion 
in molecules. In other words, interaction of DM 
clump with atomic and molecular systems is in the 
adiabatic regime. 
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